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ABSTRACT: We report the first observation of the charge symmetry breaking d + d — ► 
4 He + 7T° reaction near threshold at the Indiana University Cyclotron Facility. Kinematic 
reconstruction permitted the separation of 4 He + 7r° events from double radiative capture 
4 He + 7 + 7 events. We measured total cross sections for neutral pion production of 
12.7 + 2.2 pb at 228.5 MeV and 15.1+3.1 pb at 231.8 MeV. The uncertainty is dominated 
by statistical errors. 



Charge symmetry is the approximate symmetry of the strong interaction under a 
specific isospin rotation that interchanges down and up quarks [1,2]. This symmetry is 
broken by the different masses of the down and up quarks (md > m u ) and by their 
electromagnetic interactions. The combination of these two mechanisms leads, for example, 
to the neutron being heavier than the proton. Within the framework of chiral effective 
field theory [3,4], additional experimental information on the relative contributions of these 
two mechanisms to charge symmetry breaking (CSB) must, in leading order, come from 
pion-nucleon scattering. Direct experimental evidence is restricted to elastic scattering 
and charge exchange experiments with low-energy charged pions where the interpretation 
is complicated by corrections for the neutron-proton mass difference and electromagnetic 
interactions between the pions and nucleons [5,6]. Reactions in which a 7r° is emitted after 
being created by one nucleon and rescattered by a second are particularly clean. One 
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example of such a CSB process is the measurement of a forward-backward asymmetry in 
the cross section for the n + p — > d + n° reaction [7] . 

The d + d — > 4 He + 7r° reaction violates charge symmetry because the 7r°, whose 
wavefunction is odd under the interchange of down and up quarks, should not be produced 
from charge symmetry even (in this case self-conjugate) nuclear systems. The amplitude 
for CSB is weaker than a similar charge symmetry conserving amplitude by about 1/300, 
roughly the ratio of the quark mass difference to the nucleon mass. This suggests a 
d + d — > 4 He + 7T° total cross section, which depends only on a CSB amplitude, that is as 
small as tens of picobarns. 

Several searches for the d + d — > 4 He + ir° reaction have produced only upper limits 
[8]. A positive report at a deuteron energy of 1.1 GeV [9] has been questioned because 
the experiment did not clearly distinguish the photons from tv° decay from photons that 
could have been produced by the double radiative capture d + d — > 4 He + 7 + 7 process. 
This process is isospin allowed and predicted to appear at the reported cross section [10]. 
A further search that can clearly distinguish between these two reactions is therefore 
warranted. 

We chose to look for the d + d — > 4 He + n° reaction just above its threshold at 
225.5 MeV to avoid other pion producing channels and to take advantage of the clean 
experimental conditions afforded by the Indiana University Cyclotron Facility's electron- 
cooled storage ring. A 6° bend located in one section of the ring provided a site where 
4 He nuclei, produced in a narrow forward cone just above threshold, could be separated 
from the circulating deuteron beam. By placing a gas jet target sufficiently upstream of 
the 6° magnet, it became possible to cover a large solid angle with two arrays of Pb-glass 
detectors that would be selectively sensitive to photons from the target region. 

The layout of the experiment is shown in Fig. 1, with the major features being two 
arrays of Pb-glass detectors stacked on the left and right sides of the gas target box, the 
6° separation magnet, and a magnetic channel consisting of a septum magnet to steer the 
4 He nuclei away from the downstream ring quadrupole magnets and a quadrupole triplet 
to confine the 4 He nuclei within the acceptance of the detectors at the end of a long, 
evacuated drift length. 

For some 4 He + 7r° events, the 4 He nucleus and both 7r° decay photons could be 
recorded. Parts of the solid angle above and below the target were required for differential 
pumping of the jet target volume. Kinematic constraints just above threshold require very 
similar photon opening angles and energies for 4 He + 7r° and 4 He + 7 + 7 events. Separation 
of these two reactions relied on measuring the 4 He momentum vector using time of flight 
in the channel for the longitudinal component and scattering angle for the transverse 
component. The flight time was measured from the first plastic scintillation detector 
(AE1) to the second (AE2) at the end of the channel 5.73 m away. A third scintillator (E) 
in which the 4 He nuclei stopped was required for an event trigger along with the absence 
of any signal in two additional scintillators (VI and V2). Pulse- height correlations among 
the trigger scintillators cleanly separated 4 He events. Pb-glass information was not used 
in the trigger. The scattering angle of the 4 He was recorded by a multiwire proportional 
chamber (WC1). Two additional chambers (WC2 and WC3) tracked each 4 He nucleus 
through the channel. From the 4 He momentum and the assumption of a two-body final 
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deuterium 
jet target 

Figure 1: Layout of the experimental setup showing the target, approximate locations of the 
Pb-glass arrays, and the magnetic channel in relation to a segment of the electron- cooled 
storage ring. Quadrupole magnets (Ql, Q2, and Q3), wire chambers (WC1, WC2, and 
WC3), and scintillation trigger (AE1, AE2, and E) and veto (VI and V2) detectors are 
shown. The luminosity detectors are small and consequently are omitted here. 

state, it is possible to calculate the missing mass. Double radiative capture events produce 
a broad distribution of mass values up to a kinematic limit that depends on the beam 
energy rather than a peak at the 7r° mass. 

The isospin-allowed p ± d — > 3 He + ir° reaction was used to commission the detector 
system. The calculation of 3 He (and later 4 He) momentum from channel time of flight used 
a model to describe the energy loss (scaled from Janni [11]) of the 3 He passing through the 
channel detectors and vacuum windows. The time offsets for each photomultiplier asso- 
ciated with the AE1 and AE2 detectors were adjusted empirically during the subsequent 
data analysis. For 3 He energies close to threshold and runs lasting several hours, a ir° mass 
resolution as small as FWHM = 240 keV was obtained. 

For p + d — > 3 He + 7r° events in which both tt° photons were recorded, photon energies 
were taken to be the sum of Pb-glass energies starting with the detector recording the 
highest energy deposition and including its eight nearest neighbors. With summed energies 
above a threshold chosen to remove most random events and with Pb-glass timing in the 
correct range relative to the channel trigger, the measured efficiency for detecting the two 
7T° decay photons was 0.360 ±0.001. Simulations using the GEANT Monte-Carlo program 
library [12] and the known 7r° angular distribution [13] reproduced the Pb-glass spectral 
energy shape and agreed with the measured efficiency to ±0.01 This simulation was used 
to obtain the changes in Pb-glass efficiency between the p ± d — > 3 He ± n° commissioning 
run and production running for d ± d — > 4 He ± 7r° due to alterations of the kinematics and 
center-of-mass angular distribution (assuming the d ± d — > 4 He + n° cross section to be 



3 



isotropic). Gains for the Pb-glass detectors were calibrated using cosmic ray muons. This 
calibration operated continuously during commissioning and production running. 

The size of the 3 He opening angle cone for p + d — > 3 He + n° calibrated the incident 
proton energy. Several determinations at slightly different energies were combined by 
using the frequency of the RF voltage that maintained the beam bunching to calculate the 
circumference of the storage ring, yielding a value of 86.786 ± 0.003 m that was consistent 
over time. This calibration was used subsequently for the deuteron beam energy. 

The deuterium target was made by directing deuterium gas from a glass nozzle cooled 
to about 40 K across the circulating storage ring beam (~ 2 mm wide). The gas flow was 
adjusted so that the beam lifetime was close to 100 s with the RF voltage and electron 
cooling on, a value that gave the largest data production rate. The circulating beam 
current was as large as 2 mA. 

The luminosity (product of beam flux and intercepted target thickness) was deter- 
mined through a separate calibrated monitoring system. For this we used d + d elastic 
scattering in the vicinity of Q c . m . = 90° by placing two scintillator detectors at 44° on 
either side of the beam. The product of d + d differential cross section and solid angle 
for this system was measured in a separate experiment using a molecular HD target. Ad- 
ditional scintillators were added to observe d + p elastic scattering. Using known d + p 
cross sections and the equality of d + p and d + d luminosities for an HD target, we ob- 
tained the d + d calibration. Cross sections for d + p elastic scattering in our energy range 
have recently become available [14] with an absolute normalization error of 5.5%. The 
scintillation detectors provided particle identification information to separate breakup and 
other backgrounds. An additional position-sensitive silicon detector used to record recoil 
nuclei from small-angle scattering provided a measurement of the jet target profile as it 
intercepted the beam. This information was used to make corrections to the scintillator 
acceptance geometry for both the d + p and d + d elastic scattering processes, and to 
track changes between calibration and production running. The average luminosity during 
production was 2.9 x 10 31 /cm 2 /s. 

Candidate d + d — > 4 He + tv° events were required to have the correct pulse height 
in the three trigger scintillators (AE1, AE2, and E), usable wire chamber information, 
and photon signals in both the left and right Pb-glass arrays with an energy sum above 
threshold and timing coincident with 4 He events in the channel. For each candidate event 
the missing mass was calculated. The results were accumulated into a spectrum, as shown 
in Fig. 2. Investigations of the quality of the data set showed that random background 
was essentially removed (< 1 event /spectrum) ; thus we should interpret these results as 
arising only from 7T° production and double radiative capture. 

Production running started at 228.5 MeV, an energy chosen because the half angle 
of the 4 He forward cone (1.2°) was expected to fit safely inside the channel acceptance. 
Over the much longer d + d — > 4 He + tv° production running times, drifts in the time-of- 
flight measurement became more of a problem, despite efforts to track timing changes by 
monitoring other particles going along the channel. The width of the missing mass peak 
was FWHM = 510 keV, as seen in the top panel of Fig. 2. This meant that the peak 
was not well separated from the 4 He + 7 + 7 kinematic upper limit of 136.4 MeV. Below 
the peak, the flux we would attribute to double radiative capture is attenuated by the 
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Figure 2: Historgrams of the can- 
didate events at the two deuteron 
bombarding energies as a function 
of their missing mass value. The 
smooth curves show the reproduc- 
tion of these histograms with a 
Gaussian peak and a continuum. 
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acceptance boundaries of the channel, increasing the ambiguity of its shape. The decision 
was made to complete production running at the higher energy of 231.8 MeV (1.75° cone) 
with a kinematic endpoint of 138.0 MeV even though an additional 10% of the 7r° events 
would be lost in the channel. The result shown in the lower panel of Fig. 2 still has a 
broad peak (FWHM = 660 keV) but there is now a clear distribution of double radiative 
capture events on either side. 

The spectra of Fig. 2 were modelled with a Gaussian peak and a continuum. The 
continuum shape was obtained from the distribution in missing mass of all 4 He events (rate 
about 10 3 higher than for candidate events alone) reduced by the ratio of the calculated 
double radiative capture cross section to the phase space value. These 4 He events without 
coincident photons, which may have originated from (d, 4 He) reactions on residual gas and 
storage ring structures, were broadly distributed in energy and angle, and were used as an 
estimate of the product of phase space and channel acceptance. The continuum curves in 
Fig. 2 are a smooth representation of this shape. The centroids of the Gaussian peaks have 
a fitting error of less than 60 keV and are consistent with the tt° mass. Events in the peak 
appear to be isotropically distributed in the center of mass as shown by their distribution 
in scattering angle and time of flight. 

The 66 and 50 4 He + n° events recorded at the two energies of 228.5 and 231.8 MeV 
lead to total cross section values of 12.7 ±2.2 pb and 15.1 ±3.1 pb respectively, including a 
6.6% normalization error for all systematic effects. Corrections at the two energies included 
Pb-glass efficiency (0.34 and 0.32), trigger losses (0.94 and 0.96 for random vetoes), system 
livetime (0.95 and 0.94), wire chamber efficiency (0.93 and 0.95), and other channel losses 
from acceptance and multiple scattering (0.95 and 0.81). These cross sections are consistent 
with being proportional to r\ = Pir/m^ with a combined slope of o"tot/?7 = 80 ± 11 pb. 
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The integral of the double radiative capture process for the 2 MeV range just below the 
kinematic limit is 6.9±0.9 pb and 9.5±1.4 pb, values that are more than double predictions 
[15]. 

These results provide the first unambiguous measurement of the d + d — > 4 He + 7r° 
cross section as well as that for the competing double radiative capture process. A detailed 
interpretation of our results will require a careful consideration of different CSB reaction 
mechanisms. On the basis of chiral power counting arguments, the leading terms are 
expected to be [16] those related to n — r\ mixing [17] and CSB n + N scattering [3]. 
Preliminary plane wave calculations [16] suggest that tc — r\ mixing may be important 
because all four nucleons contribute coherently, especially if the n is produced by two 
nucleons via the exchange of a heavy meson, as in p + p — » p + p + tv° [18]. The effects 
of pion scattering are also expected to be enhanced by the initial-state d + d interactions 
[16]. Preliminary estimates show that electromagnetic effects are very small. 

Both 7T — rj mixing and CSB tv + N scattering are also important in producing an 
asymmetry about 90° in the n + p — > d + n° reaction [4] , with the latter determining the 
sign of the asymmetry. Calculations now underway should allow the extraction of both 
processes from the combined experimental results. Since both depend on the quark mass 
difference, experimental information that separates each contribution will greatly enhance 
efforts to relate the hadronic matrix elements to the underlying QCD theory. 
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